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Abstract. Our work addresses the problem of generating maximally entangled 
two spin-1/2 (qubit) symmetric states using NMR, NQR, Lipkin-Meshkov-Glick 
Hamiltonians. Time evolution of such Hamiltonians provides various logic gates which 
can be used for quantum processing tasks. Pairs of spin-l/2's have modeled a wide 
range of problems in physics. Here we are interested in two spin-1/2 symmetric states 
which belong to a subspace spanned by the angular momentum basis {\j = l,/i);/i = 
+1,0,-1}. Our technique relies on the decomposition of a Hamiltonian in terms of 

£^ ' SU(3) generators. In this context, we define a set of linearly independent, traceless, 

Hermitian operators which provides an alternate set of SU(n) generators. These 

rji matrices are constructed out of angular momentum operators J x ,J y ,J z - We construct 

and study the properties of perfect entanglers acting on a symmetric subspace i.e., 
spin-1 operators that can generate maximally entangled states from some suitably 

£> ■ chosen initial separable states in terms of their entangling power. 
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1. Introduction 

In the last few years there has been considerable increase in experimental activity [1] 
aiming to create entangled quantum states which have potential applications in quantum 
information processing tasks. In practice, these states are created by some physical 
operations involving the interaction between several systems. Thus analyzing these 
operations with regard to the possibility of creating maximally entangled states from 
an initial unentangled one and characterization of entangling capabilities of quantum 
operators play an important role in quantum information theory. Spin-1/2 (two level) 
systems have modeled a wide range of problems in physics. Considering N spin-l/2's (N 
qubits) in the symmetric subspace - the set of those N-particle pure states that remain 
unchanged by permutations of individual particles [2j [3], we define a set of (N + l) 2 
linearly independent, experimentally realizable cartesian tensor operators which provide 
different logic gates for quantum computation. In the particular case of two spin-l/2's 
NMR , NQR provide "hardware" for realizable quantum computers which involve the 
study of time evolution of Hamiltonian which can also be time dependent. Since these 
two qubit symmetric gates are capable of producing entanglement, quantifying their 
entangling capability is very important. Makhlin [I] has analyzed nonlocal properties of 
two-qubit gates and also studied some basic properties of perfect entanglers which are 
defined as the unitary operators that can generate maximally entangled states from some 
suitably chosen separable states. Zanardi et al. [5] have explored the entangling power 
of quantum evolutions in terms of mean linear entropy produced when unitary operator 
acts on a given distribution of pure product states. Kraus and Cirac |6] , Rezakhani [7] 
have given the tools to find the best separable two qubit input orthonormal product 
states such that some given unitary transformation can create maximally entangled 
quantum states. The entangling capability of a unitary quantum gate can be quantified 
by its entangling power e p (U) [5]. Balakrishnan et al. [H] have derived e p (U) in terms 
of local invariant G\. In this paper, we show that the two qubit symmetric quantum 
gates expressed in terms of newly defined linearly independent cartesian tensor operators 
belong to the class of perfect entanglers which can generate maximally entangled states 
from some suitably chosen product states. Further we show that these symmetric two 
qubit gates belong to a family of special perfect entanglers under certain conditions. 
This is a very relevant problem not only from the theoretical point of view, but also 
from the experimental one. 

In section 2, we define an alternate representation of SU(n) generators using the well 
known spherical tensor operators. Explicit form of these generators in 3-dimensional 
representation is given in 2.1. Section 2.2 deals with algebra of these cartesian tensor 
operators. Section 3 deals with two qubit symmetric gates and their entangling power 
in terms of the invariant G%. In section 3.1, we identify the conditions under which the 
perfect entangler can be classified as special perfect entangler. Further the entangling 
property of Lipkin-Meshkov-Glick Hamiltonian is studied in the spin-1 subspace. 
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1.1. Symmetric states 

Symmetric states offer elegant mathematical analysis as the dimension of the Hilbert 
space reduces drastically from 2 N to (N + 1), when N qubits respect exchange 
symmetry. Such a Hilbert space is considered to be spanned by the eigen states 
{\j, m); —j>m< +j} of angular momentum operators J 2 and J z , where j : = y. The 
corresponding density matrix gets transformed to a 3 x 3 block form in the symmetric 
subspace charactrized by the maximal value of total angular momentum j max = 1. The 
symmetric subspace provides a convenient, albeit idealized, computationally accessible 
class of spin states relevant to many experimental situations such as spin squeezing. 
Completely symmetric systems are experimentally interesting, largely because it is often 
easier to nonselectively address an entire ensemble of particles rather than individually 
address each member. Permutationally symmetric states are useful in a variety of 
quantum information processing tasks and a class of these states have recently been 
implemented experimentally [9l ITU]. 

2. Alternative representation of SU(n) generators 



^ = 727 + iy £ E h k q rf(f), (1) 



It is well known that any Hermitian operator for a spin j system is given by [H] 

1 2j +k 

k=0 q=-k 

where r q 's (with Tq = I , the identity operator) are irreducible spherical tensor 
operators of rank 'k' in the 2j+l dimension spin space with projection 'q' along the 
axis of quantization in the real 3-dimensional space. The r q 's satisfy the orthogonality 
relation 

Tr(T?T$) = (2j + l)8 kk ,6 qq ,. (2) 

Here the normalization has been chosen so as to be in agreement with Madison 
convention [12]. The spherical tensor parameters h k which characterize the given 
Hermitian operator TL are given by hJl = Tr(9iTg). Since TC is Hermitian and 
r^ = (— l) Q T^ q , hq's satisfy the condition h** = (— l) q h k _ q . The spherical tensor 
parameters h^'s have simple transformation properties under co-ordinate rotation [13] 
in the 3-dimensional space, i.e., 

(K) R = E D k q>q {aM) h k , , (3) 

q =—k 

where D k , (a/37) denote Wigner-D matrix parametrized by Euler angles (a/37). 

Following the well known Weyl construction [13] for r k in terms of angular 
momentum operators J x , J y and J z , we have 

rZ(J)=M kj (J-i) k r k Y q k (r), (4) 



Entangling capabilities of symmetric two qubit gates 

where 

2 k 



Nkn = TTi 



An(2j - k)\(2j + I) 



Jkj ~ k\\ (2j + k + l)\ 

are the normalization factors and Y k (r) are the spherical harmonics. Under rotations 
Tg's transform according to Wigner-D matrices, i.e., 

2j +k 

{r k q {J)) R = £ £ D k , q {aM) r$(J) . (6) 

fc=0 q=-k 

We now define a set of linearly independent, traceless (except (T°)q), orthonormal 
Hermitian basis matrices (T a ) k , where a = +, — ,0 , k = 1...2J , and q = 1 to +k as 
follows: 

^k _l {^k\\ 

(7) 
(8) 



(T + ) k q - 
(T~) k = 


. r k q + (r*)t 


\/2(2j + 1) ' 
_ *(?? - (r g fc )t) 



2(2j + 1) 



and 



Observe that these matrices satisfy the relation Tr{{T a )^(T^) k ,) = &ap&kk'Kq- ^ n our 
new representation the most general density matrix can be written as 

p = (r°)8(T°)g + E (^°)o(T )o fe + E E E (OjcnJ (io) 

k=l....2j a=+- k=1...2j q=l...k 

Apart from (T°)q which is proportional to identity matrix, there are 2j diagonal matrices 
namely {T°)q , k = 1...2J and the rest are off diagonal. 

2.1. SU(3) generators: 

In the particular case of two qubit symmetric subspace, our set of basis matrices |§| can 
be obtained from Equation(7,8,9) asljj] 

Mo-^ , M 1 = — ? =-, M 2 = ^ , 

M 3 = W-r 1 , M 4 = l 2 _ 2 ; , M 6 - ' 



3 U ' "*• v^ ' ° _V / 3 ' 

r 2 4- r 2t r 2 4- r 2t /2 

M 6 = ^, M 7 ^^, M 8 ^. (H) 

§ The above matrices are equivalent to the set of matrices with different normalization defined by 

R.J.Morris [15]. 

|| We have used a different notation for the basis set for the sake of simplicity. 
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These operators are explicitly represented in |lm) basis where m 
follows: 
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The above matrices are normalized i.e., Tr(M k M k i) = 2 5 kk > and Mi,...,M 7 have 
eigen values 1, 0, -1. For a pair of spin 1/2's (two qubits), there are 4x4 linearly 
independent operators that close under mutual commutator brackets. The 16 linearly 
independent operators of a four-state system can be chosen in a variety of matrix 
representations. One such choice [20] is |, |, and §<E> | where a and r are the individual 
pauli matrices for two qubits. Together with the 4x4 unit matrix, these 16 operators can 
be used to construct any operator describing magnetic couplings between the spins as 
well as the coupling of each spin to an external field. The co-efficient of these operators 
may in general be functions of time and have physical significance. In NMR and quantum 
computation, the effect of such operator combinations on paired spins and the solutions 
of the corresponding Hamiltonians is considered. Our complete set of basis matrices 
which are 9 in number are unitarily equivalent to combinations of these operators, i.e., 

P 0l — ► M , -(O s + 9 ) — ► M 1 , (O e + do) 



(02 + 03) — 

-2(0 7 + O n ] 



M 3 , -2(Ois + Oie) — ■> M 4 , 
— ->M 6 ,2(Oi3-Oi 4 )-^M 7 



^M 2 

2(08 + Oi 2 ) 
2 



Ms 



V3 



(20 4 - 13 - 



uj 



M 



8 • 



And the unitary transformation which takes an operator from qubit basis to angular 
momentum basis is 

0\ 



u 



( 1 











1 

v/2 


1 

V2 











V° 


1 


1 

V2 



As |lm) basis is related to the qubit basis through 1 11) = |^f) , 1 10} 
|1 — 1) = \l\) , the above 9 matrices in the qubit basis are realized as 



^±p,and 



Mo = t/|(|tt)(ttl + \U)(U\) + J((in> + Utl) + ((til +(ltl)) , 
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Mi = - i(|tt)((tll +UTD+ ltl)((ttl +(111) 
+ llt)((ttl +(U\)+ \UK(U\ +UTD) , 

m 2 = ^(|tt)((tll +(ITI)+ ItDHTTI +(111) 
+ IITX-Crtl +(U\)- \UK(U\ +(ITD) , 

m 3 =(|tt)(ttl)-(lll)(lll), 

m 4 =z((|||}(tt|)-(|tt)ail)), 

m 5 = Jattxmi +ati)- iti)((tti +aii) 
-iit)((tti+aii)+iii)((tn+ati)), 

m 6 = l -{- itt)((tn +ati)+ itix-crti +uii) 
+ iit)(-(tti+(ui)+ni)((tii+ati)), 

M 7 = ((itt)aii) + (iii)(ttD), 

m 8 = -L(dttxtti)- wxmi +ati)- nt)((tn +ati)+ iid(ui) • 

In this representation the most general spin-1 Hamiltonian can be written as 

1 8 
n(t) = -J2 h k (t)M k . (12) 

1 i=0 

Here M^'s in terms of angular momentum operators J x ,J y ,J z are given by Mi = 
-(J x ) , M 2 = (J y ), M 3 = (J 2 ), M 4 = -(J x J y + J y J x ) , M 5 = (J y J z + J 2 J y ) , 
M 6 = -(J X J 2 + J 2 J x ) , M 7 = (Jl - J^) , M 8 = (3J 2 2 - 2) . Note that the expansion 
co-efficients hk = Tri^-iM^) are real and hence they constitute an experimentally 
measurable set of parameters. 

2.2. Algebra of Cartesian tensor Operators 

Commutation and anticommutation relations between M^'s are given in table I and II 

respectively. 

Table I: Table of commutators [Mj,, M#], with operators M& in the first column and M^ 
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in the top row, each entry provides the commutator [M k , M k >]. 



M k 
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VSiM 5 


M 2 


iM 3 
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M 4 


iM 6 


-iM 5 


-2iM 7 
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-vL, (13 > 
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2iM 4 


-2iM 3 
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VSiM 6 








v^^Mi 


-V3iM 2 









where a = i(>/3Af 8 + M 7 ) , b = i(VSM 8 - M 7 ). Here [M 3 ,M 8 ]=0, [M 4 ,M 8 ]=0, 
[M 7 ,M 8 ]=0. Also there are four vector triplets [Mi,M 2 ,M 3 ], [Mi,M 4 ,M 6 ], 
[M 4 ,M 2 ,M 5 ], [M 5 ,M 3 ,M 6 ] satisfying the relation [Mi.Mj] = -ie ijk M k and one triplet 
[M 4 ,M 3 ,M 7 ], satisfying -2ie ijk M k . 

Table II: Table of anticommutators {M k , M k /}, with operators M k in the first 
column and My in the top row, each entry provides the anticommutator {M k ,M k >}. 
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"73 Ml 


-73 M2 
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3s M * 


"75^ 


-73 M « 


T^ 


L> 



(14) 



A = 2,/|m + M 7 - 4=M 8 , B = 2,/|m - M 7 - -J-M 8 , C = 2,/|m + 4.M f 



3^0 -r,„ 7 ^« 8 ,^-^ V3 ,™ .« 7 ^ H8 ,v- iV 3^ T ^ 



D = 2 A /|Mn-^M«. 



3 J "° - 7^8 ■ 
The relationship between M k 's, and the Gell-Mann matrices A' fe s k=1...8 is given 

by 

Mi = "-4(Ai + A 6 ), M 2 = -J=(A 2 + A 7 ), M 3 = ^A 3 + ^A 8 , M 4 = -A 5 , 

M 5 = -^=(A 2 - A 7 ), M 6 = -^=(A 6 - Ai), M 7 = A 4 , M 8 = ^A 3 - ^A 8 . 
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3. Two qubit symmetric gates 



The algebraic study of M' k s and the corresponding symmetric two qubit gates is 
important not only for understanding fundamental properties of coupled spins and of 
quantum circuits, but also to study possible experimental implementations in different 
physical systems. Hamiltonian evolution provides the hardware for quantum gates. 
i.e., the time evolution of the operators M^'s provide various symmetric logic gates for 
quantum computation. The closed form expression for e lMk9 is given by Bk = e tMke = 
I + (cos9 — 1)M| + isinOMk- Here k = 1....7 and I is a 3 x 3 unit matrix. Following are 
the explicit forms of the gates B^s in the symmetric subspace: 
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A useful property of a two qubit symmetric gate is its ability to produce a maximally 
entangled state from an unentangled one. This property is locally invariant. It is well 
known that perfect entanglers are those unitary operators that can generate maximally 
entangled states from some suitably chosen separable states. The entangling properties 
of quantum operators have already been discussed in the literature [5J [8j [H] . Here we 
calculate the entangling power of two qubit symmetric gates following the simplified 
expression given by Balakrishnan et al. [8] according to which the gate B is a perfect 
entangler if its entangling power, e p (B) = |(1 — |Gi|) has the range | < e p < |. 
The local invariant G\ [Ref. [1] table II] in terms of symmetric, unitary matrix m is 

Here m = B^Bb where the gates in the Bell basis are given by 



tr^m 



given by G x - wdet[B] 

UBU^ . U is a transformation matrix given by 



B 



B 






( 1 





1 





-V2i 














V -i 
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\ 



1), 1 00) to the Bell basis lu> ^ n> . 



connecting the angular momentum basis |11), |10), 

K1M£M 1-ltHt-l) ^(l-UHtt)) _ The relation 6p ^ = |(i _ | dl) implies that gates having 

the same |Gi| must necessarily possess the same entangling power e p . 

It is obvious that B\, B2, B3 donot produce entanglement as they represent rotations 
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which is a local unitary transformation. Note that |C?i| = 1 and e p = for the above 
gates. Interestingly, for the gates S 4 , B 5 , B 6 and B 7 , \Gi\ = Cos A {9). Observe that 
since < G x < 1 for < < § , it is clear that < e p (B B ) k < § ( k = 4.. .7). All these 
above mentioned gates are perfect entanglers for ~ < 6 < ^. Similarly the gate Bg will 
have maximum entangling power i.e., e p = 2/9 when 9 = y3| 



As an example, consider the direct product state \1p12) 
in the qubit basis. 



\i/ji) <S> \ip2) , of two spinors 
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/ 
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< «i,2 < 7T , < 0i : 2 < 27T . Note that a separable state in the symmetric subspace 
will have the form 
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\ 
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where «i = a 2 = a and 4>\ = 2 = 4>- 

It is a well known fact that for a pure state of two qubits \ tfj) = a |tt) + b |tl) + c lit 
) + d \ll), the expression for concurrence is C(ip) = 2\ad — bc\ [15J. For a maximally 
entangled quantum state concurrence C = 1. It can be observed that under the action 
of the gates i? 4 , B 7 and B 8 (with e p being maximum i.e., 2/9), \ipi2) S ym will become 
maximally entangled state when a = ~. i.e., 
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or in the qubit basis 
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Similarly, the gates B 5 , B e acting on the symmetric separable state transform it into 
maximally entangled one when a = 0, n. For eg: 
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It can be noted that the operators B 8 and B4 produce spin squeezing resulting from a 
single axis twisting and two axis counter twisting respectively [16]. Also possibility of 
physical realization of these spin squeezing operators are given in Ref.|18j. 



3. 1 . Special perfect entanglers 

Rezakhani [7] has analyzed the perfect entanglers and found that some of them have 
the unique property of maximally entangling a complete set of orthonormal product 
vectors. Such operators for which e p — | belong to a well known family of special 
perfect entanglers. A study of using such special perfect entanglers as the building 
blocks of the most efficient universal gate simulation is also given in ref.[7]. Let us 
now study the conditions under which the perfect entanglers can be classified as special 
perfect entanglers. When e p = |, _B 4 , B 8 in the qubit basis are given by 
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Following Rezakhani [7], the most general separable basis (upto general phase factors 
for each vector) is 



IV>i> 
IV*> 

1^3) 
1^4) 

where \a\ 2 + \b\ 2 = 



(a\t) + b\l))®(c\t} + d\l)), 

(-6'|t>+a*||»®(c|t> + dU», 

(e|t) + /U))®M1t) + c1l)), 
(-/*|t) + e*U))®M*|t)+c*U)), 



+ \d? 



+ I/P 



1. 
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When the gates B±, B 7 and B s as perfect entanglers act on the state - say \ipi), we 
obtain 

[#4,7,8] 1-01 ) = -bd | ft) + ad | t!) + be | ||> + ac | ||). 

This state is maximally entangled if its concurrence, C = 4|a6c<i| = 1. Thus these two 
qubit symmetric gates transform the orthonormal states \tl>i), \i>2), I ^3) an d l^) m t° 
maximally entangled ones if \abcd\ = \cdef\ = \. Similarly, for the gates B 5 and B 6 , 
condition for finding a full set of orthonormal product states is \(a 2 + b 2 )(c 2 + d 2 )\ = 
\(e 2 + f 2 )(c 2 + d 2 )\ = 1. 

Let us consider the example of Lipkin-Meshkov-Glick interaction Hamiltonian 
[T7J [18] which is widely used in nuclear physics. 

n L = g 1 (j 2 + + j 2 _) + g 2 (J + J- + J-J + ) . (15) 

Here Q\ and Q2 are the coupling constants. In terms of our operators Mk 's, 

n L = g[M 7 + g' 2 (V8M - m 8 ) , (16) 

where Q[ = 2Qi and Q' 2 = -jnQ2- Since [M 7 , M 8 ] = 0, we have 

/ e V3iP C0S {: ie^sin^ \ 
e lH ^ = B L = e 2 ^ 

v ie^ ifi cosi e^^cosi j 

in spin-1 subspace. Here ^ = Q[t and j3 = Q' 2 t and e p — I for 2Q 2 t = f + 2Qit. Under 
the action of this gate (with e p = |), the separable state | ft) (I ID) becomes entangled 

for all values oft except when t = j£-; n=0,l,2 and maximally entangled when AQit 

= (2n + l)f. For eg., 

a = 
B L \Msvm — > cos{2Q 1 t) Ift) +iain(2g l t) \U) ■ 

4. Conclusion 

In conclusion, we have constructed traceless,Hermitian and linearly independent set 
of basis matrices which provides an alternative representation of SU(n) generators. We 
have considered unitary evolutions of two spin-1/2 states in angular momentum subspace 
(j=l) and constructed physically realizable logic gates using (2j+l) dimensional 
representation of the above set of basis matrices. Entangling properties of these gates 
have been studied in terms of their entangling power e p . e p is found to be maximum 
(2/9) for B^,...,B S under certain conditions which is the signature for special perfect 
entanglers. These logic gates are obtained by the exponentiation of the quadratic form 
of angular momentum operators J x ,J y ,J z . As an example we have taken the well 
known Lipkin-Meshkov-Glick Hamiltonian and studied its entangling properties in spin- 
1 subspace. Further, we have shown that precisely at what time the initial separable 
state becomes maximally entangled under the action of perfect entanglers which consists 
of one-axis twisting and two axis twisting Hamiltonians that produce spin squeezing. 
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